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Abstract — We consider distributed consensus and vehicular 
formation control problems. Specifically we address the question 
of whether local feedback is sufficient to maintain coherence in 
large-scale networks subject to stochastic disturbances. We define 
macroscopic performance measures which are global quantities 
that capture the notion of coherence; a notion of global order 
that quantifies how closely the formation resembles a solid 
object. We consider how these measures scale asymptotically 
with network size in the topologies of regular lattices in 1, 2 
and higher dimensions, with vehicular platoons corresponding to 
the 1 dimensional case. A common phenomenon appears where 
a higher spatial dimension implies a more favorable scaling of 
coherence measures, with a dimensions of 3 being necessary to 
achieve coherence in consensus and vehicular formations under 
certain conditions. In particular, we show that it is impossible to 
have large coherent one dimensional vehicular platoons with only 
local feedback. We analyze these effects in terms of the underlying 
energetic modes of motion, showing that they take the form of 
large temporal and spatial scales resulting in an accordion-like 
motion of formations. A conclusion can be drawn that in low 
spatial dimensions, local feedback is unable to regulate large- 
scale disturbances, but it can in higher spatial dimensions. This 
phenomenon is distinct from, and unrelated to string instability 
issues which are commonly encountered in control problems for 
automated highways. 



I. Introduction 

The control problem for strings of vehicles (the so-called 
platooning problem) has been extensively studied in the last 
two decades, with original problem formulations and studies 
dating back to the 60's |[T|-||5|. These problems are also 
intimately related to more recent formation flying and for- 
mation control problems ||6J. It has long been observed in 
platooning problems that to achieve reasonable performance, 
certain global information such as leader's position or state 
need to be broadcast to the entire formation. A precise 
analysis of the limits of performance associated with localized 
versus global control strategies does not appear to exist in 
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the formation control literature. In this paper we study the 
platooning problem as the 1 dimensional version of a more 
general vehicular formations control problem on regular lat- 
tices in arbitrary spatial dimensions. For such problems, we 
investigate the limits of performance of any local feedback 
law that is globally stabilizing. In particular, we propose and 
study measures of the coherence of the formation. These are 
measures that capture the notion of how well the formation 
resembles a rigid lattice or a solid object. 

The coherence of a formation is a different concept from, 
and often unrelated to, string instability. In the platooning 
case (i.e. 1 dimensional formations), which turns out to be 
most problematic, a localized feedback control law may posses 
string stability in the sense that the effects of any injected 
disturbance do not grow with spatial location. However, as we 
show in this paper, it is impossible to achieve a large coherent 
formation with only localized feedback if all vehicles are 
subject to any amount of distributed stochastic disturbances. 
The net effect is that with the best localized feedback, a 1 
dimensional formation will appear to behave well on a "micro- 
scopic" scale in the sense that distances between neighboring 
vehicles will be well regulated. However, if a large formation 
is observed in its entirety, it will appear to have temporally 
slow, long spatial wavelength modes that are unregulated, 
resembling an "accordion" type of motion. This is not a safety 
issue, since the formation is microscopically well regulated, 
but it might effect throughput performance in a platooning 
arrangement since throughput does depend on the coherence 
or rigidity of the formation. 

The phenomenon that we discuss occurs in both consensus 
algorithms and vehicular formation problems. We therefore 
treat both as instances of networked dynamical systems with 
first order and second order local dynamics respectively. Both 
problems are set up in the d-dimensional torus Z^. We begin 
in section |ll] with problem formulations of the consensus type 
and vehicular formations, where we view the former as a 
first order dynamics version of the latter In section [III] we 
define macroscopic and microscopic measures of performance 
in terms of variances of various quantities across the net- 
work. We argue that the macroscopic measures capture the 
notion of coherence. We also present compact formulae for 
calculating those measures as norms of systems with 
suitably defined output signals. These norms are calculated 
using traces of system Grammians, which in turn are related 
to sums involving eigenvalues of the underlying system and 
feedback gains matrices. Since the network topologies we con- 
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sider are built over Tori networks, these system matrices are 
multi-dimensional circulant operators, and their eigenvalues 
are calculated as the values of the Fourier symbols of the 
underlying feedback operators, thus allowing for a rather direct 
relation between the structure of the feedback gains and the 
system's norms. Much of the remainder of the paper is devoted 
to establishing asymptotic (in network size) bounds for these 
performance measures for each underlying spatial dimension. 



Section IV establishes upper bounds of standard algorithms. 



while section |V] is devoted to establishing lower bounds for 
any algorithm that satisfies a certain number of structural 
assumptions including the locality of feedback and bounded- 
ness of control effort. This shows that asymptotic limits of 
performance are determined by the network structure rather 
than the selection of parameters of the feedback algorithm. We 
pay particular attention to the role of control effort as our lower 
bounds are established for control laws that have bounded 
control effort in a stochastic sense. Some numerical examples 
illustrating the lack-of-coherence phenomenon are presented in 
section |VI| as well as an illustration of how it is distinct from 
string instability. The interested reader may initially skim this 
section which numerically illustrates the basic phenomenon 
we study analytically in the remainder of the paper We end 
in section IVIII with a discussion of related work in which 
various versions of this phenomenon were observed, as well 
as a discussion of some open questions. 

Notation and Preliminaries 

The networks we consider are built over the d-dimensional 
Torus Z^. The one-dimensional Torus Z^r is simply the set of 
integers {0, 1, . . . , iV— 1} with addition modulo N ( mod N), 
and is the direct product of d copies of Z^r. Functions 
defined on Z^ are called arrays, and we use multi-index 
notation for them, as in afc — a(i.j /j^) to denote individual 
entries of an array. Indices are added in the Z^ arithmetic as 
follows 



,rd) = (fci 



,kd) + {h,...,ld) 



ih + k) 



N ' 



i = 1, . . . 



where {)m is the operation mod A^. The set Z^ and the 
corresponding addition operation can be visualized as a "circu- 
lant" graph in d-dimensional space with edge nodes connected 
to nodes on corresponding opposite edge of the graph. 

The multi-dimensional Discrete Fourier Transform is used 
throughout. All states are multi-dimensional arrays which we 
define as real or complex vector-valued functions on the Torus 
Z^. The Fourier transform (Discrete Fourier Transform) of an 
array a is denoted with a. We refer to indices of spatial Fourier 
transforms as wavenumbers. Generally, we use k and I for 
spatial indices and n and m for wavenumbers. For example, an 
array a(j.j j.^) has (fci, . . . , kd) as the spatial index, while its 
Fourier transform a(„j has the index (711, . . . , Ud) as the 
wavenumber. The wavenumber is simply a spatial frequency 
variable. Some elementary properties of this Fourier transform 
are summarized in Appendix [A] 



Convolution operators arise naturally over Z^. Let a be any 
array of numbers (or matrices) over Z^, that is a : Z^ — > C 
(or C"^"). Then the operator of multi-dimensional circular 
convolution with the array a is defined as follows 

9 = Taj = a*/ 
t 

9(ku....ka) = 51 a.(ki,...,ki)~(lu...,li) f(h.,...,u)- 

Note that / and g may be scalar or vector-valued (depending 
on whether a is scalar or matrix-valued respectively), and that 
the arithmetic for (/ci, . . . , kd) — [li, . . . ,ld) is done in Z^, 
i.e. arithmetic mod N in each index as described above. 

It is important to distinguish between an array a and the 
corresponding linear operator Ta- The Fourier transform a 
of the array a is called the Fourier symbol of the operator 
Ta- It is a standard fact that the eigenvalues of the operator 
Ta are exactly the values of the Fourier transform a, i.e. the 
values of its Fourier symbol. When a is matrix valued, then the 
eignvalues of Ta are the union of all eigenvalues of a(„i,...,nj) 
as the wavenumber (ni, . . . , n^) runs through Z^, i.e. 



<y{Ta) 



u 



(ni,...,n<j)eZ 



where a{.) refers to the spectrum of a matrix or operator (all 
finite-dimensional in our case). 

In this paper, we use the term dimension to refer exclusively 
to the spatial dimension of underlying networks. To avoid 
confusion with the notion of state dimension, we refer to the 
dimension of the state space of any dynamical system as the 
order of that dynamical system. 

The vector dimension of signals is mostly suppressed to 
keep the notation from being cumbersome. For example, the 
state of node (fci, . . . , kd) in the d-dimensional Torus is written 
as 

X(ku-.ka){t)- 

It is a scalar-valued signal for consensus problems, and vector- 
valued (in M'') signal for vehicular formation problems. 

We use to denote the transpose of a matrix M, and 
M* to denote the complex-conjugate transpose of a matrix 
M or the adjoint of an operator M. Although all operators in 
this paper are finite dimensional, we sometimes refer to them 
as operators rather than matrices since we often avoid writing 
the cumbersome explicit matrix representations (such as in the 
case of multi-dimensional convolution operators). 

II. Problem Formulation 

We formulate two types of problems, consensus and ve- 
hicular formations. The mathematical setting is analogous in 
both problems, with the main difference being that vehicular 
models have two states (position and velocity) locally at each 
site in contrast to a scalar local state in consensus problem. 
This difference leads to more severe asymptotic scalings in 
vehicular formations as will be shown in the sequel. 
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A. Consensus with random insertions/deletions 



B. Vehicular Formations 



We begin by formulating a continuous-time version of the 
consensus algorithm with additive stochastic disturbances in 
the dynamics p), |8|. As opposed to standard consensus 
algorithms without additive disturbances, nodes do not achieve 
equilibrium asymptotically but fluctuate around the equilib- 
rium, and the variance of this fluctuation is a measure of how 
well approximate consensus is achieved. This formulation can 
be used to model scenarios such as load balancing over a 
distributed file system, where the additive noise represents file 
insertion and deletion, parallel processing systems where the 
noise processes model job creation and completion, or flocking 
problems in the presence of random forcing disturbances. 

We consider a consensus algorithm over undirected tori, 
Z^, where the derivative of the scalar state at each node is 
determined as a weighted average of the differences between 
that node and all its 2d neighbors. One possible such algorithm 
is given by 

ifc = /3 [{x(^ki-i....,ka) -Xk) -\ h 

{X(ki,...,ka + 1) - Xk)] + Wk, 

= (-2d/3) Xfe + /3 + X(^ki+i,....k^) + 

■■■+X(^ku.-k^-l) + X{k,.^...M^ + l)] + Wk, (1) 

where we have used equal weights /? > for all the differ- 
ences. The process disturbance w is a mutually uncorrelated 
white stochastic process. We call this the standard consensus 
algorithm in this paper since it is essentially the same as other 



well-studied consensus algorithms ||9|-|[T3|. 

The sum in the equation above can be written as a 
multimulti-dimensional convolution by defining the array 



-2d/3 



O 



(ki,...,kd) 







fcl = • • • = 

h = ±i, 

otherwise. 



= fcrf - 0, 

and kj = for i 7^ j, 



The system ([TJ can then be written as 

X — O -k X + W, 



(2) 



(3) 



where * is circular convolution in Z 

We recall that we use the operator notation Ta x := a-kx 
to indicate the circulant operator of convolution with any array 
a. With this notation, a general spatially invariant consensus 
algorithm can be written abstractly as 



Ta X + W, 



(4) 



for any array a defined over Z^. Such algorithms can be 
regarded as a combination of open loop dynamics 



Xk 



Uk 



Wk, 



with the feedback "control" u = TaX, where the feedback 
operator array is to be suitably designed. With this point of 
view, consensus algorithms can be thought of as first order 
dynamics versions of vehicular formation problems that we 
introduce next. 



Consider N'^ identical vehicles arranged in a d-dimensional 
torus, Z^, with the double integrator dynamics 

x(ki,...,ka) = U(^ki^...,kd) + (5) 

where {ki, . . . , kd) is a multi-index with each ki E Z^, u 
is the control input and it; is a mutually uncorrelated white 
stochastic process which can be considered to model random 
forcing. In the sequel, we will also consider the consequences 
of the presence of viscous friction terms in models of the form 

i(fci,....fci) = + "(fci....,/cd) + (6) 

where /i > is the linearized drag coefficient per unit mass. 

Each position vector Xk is a d-dimensional vector with 
components Xk = [ x^ ■ ■ ■ ] . The objective is 
to have the fcth vehicle in the formation follow the desired 
trajectory Xk 



A, 



which means that all vehicles are to move with constant 
heading velocity v while maintaining their respective position 
in a Z^ grid with spacing of A in each dimension. The 
situation of different spacings in different directions can be 
similarly represented, but is not considered for notational 
simplicity. 

The deviations from desired trajectory are defined as 





"4" 










Xk ■■= vt + kA ^ 


. -^k . 






t + 


_ fed . 



Xk 



Xk 



Vk 



Xk - V. 



We assume the control input to be full state feedback and 
linear in the variables x and v (therefore affine linear in x 
and v), i.e. u = Gx + Fv, where G and F are the linear 
feedback operators. The equations of motion for the controlled 
system are thus 

d 



dt 



X 







I ' 




X 




■ " 


V 




G 


F 




V 


+ 


/ 



w. 



(7) 



We note that the above equations are written in operator form, 
i.e. by suppressing the spatial index of all the variables. 

Example: The operators G and F will have some very 
special structure depending on assumptions of the type of 
feedback and measurements available. Consider for example 
a feedback control of the fcth vehicle (in a one dimensional 
formation) of the following form 



Uk 



fe+i 



f+ (vk+i 
9o [xk - 



~Xk~ A) + g_ {xk-i - Xk- A) 
-Vk) + L{vk-i-Vk) 

Xk) + fo{vk-v), 



where the g's and /'s are design constants. The first two lines 
represent look-ahead and look-behind position and velocity 
error feedbacks respectively. We refer to such terms as relative 
feedback since they only involve measurements of differences. 
On the other hand, terms in the last line require knowledge 
of positions and velocities in an absolute coordinate system (a 
grid moving at constant velocity), and we thus refer to such 
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terms as absolute feedback. For later reference, it is instructive 
to write the feedback in the above example in terms of the state 
variables x and v as 

■Uk = g+{xk+i-xk) + g_{xk-i - Xk) + 

f+{vk+l~Vk) + f.{vk-i-Vk) + 

9o Xk + fo Vk- (8) 

C. Structural assumptions 

We now list the various assumptions that can be imposed on 
system operators and on the control feedbacks G and F. These 
are structural restrictions representing the structure of open 
loop dynamics and measurements, and the type of feedback 
control available respectively. 

(Al) Spatial Invariance. All operators are spatially in- 
variant with respect to Z^. This implies that they are 
convolution operators. For instance, the operation Gx can 
be written as the convolution (over Z^) of the array x with 



an array {G(^ku.-k^)} 



(Gx) 



fe-i Xi, 



(9) 



lei 



where the arithmetic for k — lis done in Z^. For each k, the 
array element Gk is a d x d matrix (G is then an iV* x N"^ 
operator). Note that in the absence of spatial invariance, 
each term of the sum in (|9| would need to be written as 
Gkj xi. That is, one would require a two-indexed array of 
matrices Gkj rather than a single-indexed array. 
In the example above of a one dimensional circular forma- 
tion, the array elements for position feedback are given by 
{(.9o-5+ -.9-) ,5+,0, . . . ,0,.g_}. 
(A2) Relative vs. Absolute Feedback. We use the term 
Relative Feedback when given feedback involves only dif- 
ferences between quantities. For example, in position feed- 
back, this implies that for each term of the form ax(^ki....,ka) 
in the convolution, another term of the form — Q;a:(ij ... 
occurs for some other multi-index I. This implies that the 
array G has the property 



J2 Gk ^ 0. 



(10) 



fcez 



We use the term Absolute Feedback when given operator 
does not satisfy this assumption. 

Note that in the example above, relative position feedback 
corresponds to go ~ 0, and in this case, condition (lOi is 
satisfied. 

(A3) Locality. The feedbacks use only local information 
from a neighborhood of width 2q, where q is independent 
of N. Specifically, 

G(fci,...,fc^) = 0, if ki > q, and ki < N -q 

for any i e {1, . . . , d}. (11) 

The same condition holds for F. 

(A4) Reflection Symmetry. The interactions between ve- 
hicles have mirror symmetry. This has the consequence that 
the arrays representing G and F have even symmetry, e.g. 



for each nonzero term like ceG(^ki....,ka} the array there 
is a corresponding term aG(^^ki,....-kd)- This in particular 
implies that the Fourier symbols of G and F are real valued. 
In the example above, this condition gives ~ g and 

(A5) Coordinate Decoupling. For d> 2, feedback control 
of thrust in each coordinate direction depends only on mea- 
surements of position and velocity error vector components 
in that coordinate. This is equivalent to imposing that each 
array element Gk and Fk are d x d diagonal matrices. For 
further simplicity we assume those diagonal elements to be 
equal, i.e. 

Gk = diag{gfc, . . . ,.gfe}, Fk = dmg{ f k, ... Jk}- (12) 

This in effect renders the matrix-vector convolution in (|9]) 
into d decoupled scalar convolutions. 

Assumptions (Al) through (A3) appear to be important for 
subsequent developments, while assumptions (A4) and (A5) 
are made to simplify calculations. 

III. Performance Measures 

We will consider how various performance measures scale 
with system size for the consensus and vehicle formations 
problems. Some of these measures can be quantified as steady 
state variances of outputs of linear systems driven by stochastic 
inputs, so we consider some generalities first. Consider a 
general linear system driven by zero mean white noise with 
unit covariance 

X = Ax + Bw, 
y — Hx. 

Since we are interested in cases where A is not necessarily 
Hurwitz (typically due to a single unstable mode at the origin 
representing motion of the mean), the state x may not have 
finite steady state variances. However, in all cases we consider 
here the outputs y do have finite variances, i.e. the unstable 
modes of A are not observable from y. In such cases, the 
output does have a finite steady state variance, which is 
quantified by the square of the norm of the system from 
w to y 



V := 



J2 lim E{yUt)yk{t)}., 



(13) 



where the index k ranges over all "sites" in the lattice Z^. 

We are interested in spatially invariant problems over dis- 
crete Tori. This type of invariance implies that the variances 
of all outputs are equal, i.e. E {y'^yk} is independent of k. 
Thus, if the output variance at a given site is to be computed, 
it is simply the total norm divided by the system size 

E{ylyk} - ^ E ^iy^y^} = ^' (14) 



lei 



where M is the size of the system (Al — N'^ for d-dimensional 



Tori). We refer to quantities like ( 14 1 as individual output 



variances. 
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Next, we define several different performance measures and 
give the corresponding output operators for each measure for 
both the consensus and vehicular formation problems. In the 
vehicular formation problem, we assume for simplicity that the 
output involves positions only, and thus the output equation 
has the form 



y 



[c 0] 



i.e. H = [ C ], where C is a circulant operator. A 
consensus problem with the same performance measure has 
a coiTesponding output equation of the form (with the same 
C operator) 

y = Cx. 

Performance Measures: We now list the three different 
performance measures we consider. 

(PI) Local error. This is a measure of the difference 
between neighboring nodes or vehicles. For the consensus 
problem, the fcth output (in the case of one dimension) is 
defined by 

Vk ■= Xk - Xk-i- 

For the case of vehicular formations, local error is the 
difference of neighboring vehicles positions from desired 
spacing, which can equivalently be written as 

Vk ifc - ifc-i- 

The output operator is then given by C :— {I — D), where 
D is the right shift operator, {Dx)k '■— Xk-i- 
In the case of d dimensions, we define a vector output that 
contains as components the local error in each respective 
dimension, i.e. 



C 



I - D" 



(15) 



where is the right shift along the rth dimension, i.e. 

P''a;)(fci,...,A,v,...,fcd) x^ki,...,k^-u...,ka), and 1/V2d is 
a convenient normalization factor. This operator is closely 
related to the standard consensus operator O in Eq. Q by 
the following easily established identity 



C*C 



-1 
2dl3 



O. 



(16) 



(P2) Long range deviation (Disorder). In the consensus 
problem, this corresponds to measuring the disagreement 
between the two furthest nodes in the network graph. As- 
sume for simplicity that N is even and we are in dimension 
1. Then, the most distant node from node fc is ^ hops away, 
and we define 



Vk ■= Xk 



^k+i 



In the vehicular formation problem, long range deviation 
corresponds to measuring the deviation of the distance 
between the two most distant vehicles from what it should 
be. The most distant vehicle to the fcth one is the vehicle 



indexed by fc 



N 



The desired distance between them is 



A^, and the deviation from this distance is 



Vk 



Xk 



•-k+s^ 



A- 



N 



Xk 



(17) 



We consider the variance of this quantity to be a measure 
of disorder, reflecting the lack of "end-to-end rigidity" in 
the vehicle formation. 

Generalizing this measure to d dimensions yields an output 
operator of the form 



C 



- (50_5(N/2,- ■,«/2)j, (18) 

i.e. the operator of convolution with the arra}^' S'^ — 

giN/2,....N)2)^ '-' 

(P3) Deviation from average. For the consensus problem, 
this quantity measures the deviation of each state from the 
average of all states. 



Vk ■= Xk 



M ^ 



Xi. 



(19) 



In operator form we have y = (/ — Ti)x, where 1 is the 
array of all elements equal to 1/M. 

In vehicular formations, this measure can be interpreted 
as the deviation of each vehicle's position error from the 
average of the overall position error y — {I — Ti)x. 
We note that performance measures (PI) through (P3) are 
such that C can be represented as a convolution with an array 
{Ck} which has the property ^ ^.^^d Ck = 0. This condition 
causes the mean mode at zero to be unobservable, and thus 
guarantees that all outputs defined above have finite variances. 

We refer to the performance measure (PI) as a microscopic 
error since it involves quantities local to any given site. This is 
in contrast to the measures (P2) and (P3) which involve quanti- 
ties that are far apart in the network, and we thus refer to these 
as macroscopic errors. We consider the macroscopic errors as 
measures of disorder or equivalently, lack of coherence. As 
we will show in the sequel, both macroscopic measures scale 
similarly asymptotically with system size, which justifies using 
either of them as a measure of disorder. 

Formulae for variances: Since we consider spatially invari- 
ant systems and in particular systems on the discrete Tori Z^, 
it is possible to derive formulae for the above defined measures 
in terms of the Fourier symbols of the operators K, F and C. 
Recall the state space formula for the norm V defined 
in ^ 



V = tr 



B*e^'*H*He^*B dt 



When A, B and H are circulant operators, traces can be 
rewritten in terms of their respective Fourier symbols (see 
(|47])) as 



V = 



inJ2 Ke^'-'HlH^e^-'Bn dt (20) 



(21) 



'By a slight abuse of notation, we define tlie sliifted Kronecker delta 5 J. : = 
&k-l^ where (Sj. = 1 for = 0, and zero otherwise, is the standard Kronecker 
delta. With this notation, 5" is also the standard Kronecker delta. 
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where the individual integrals are defined as 



Pn ■■ = 



-,Alt 



(22) 



If An is Hurwitz, then _P„ can be obtained by solving the 
Lyapunov equation 



KPn 



P A 



(23) 



For wavenumbers n for which An is not Hurwitz, P„ is still 
finite if the non-Hurwitz modes of A„ are not observable from 



Hn- In this case we can analyze the integral in (22 1 on a case 
by case basis. 

The Lyapunov equations are easy to solve in the Fourier 
domain. Equation ( [23] l is a scalar equation in the Consensus 
case and a 2d x 2d matrix equation in the Vehicular cas^ 
The two respective calculations are summarized in the next 
lemma. The proof is given in the Appendix. 

Lemma 3.1: The output variances 



13 I for the consensus 



and vehicular problems satisfying assumptions (A1)-(A5) are 
given by 

_ 1 ^ 1^ 12 



3fi(a„)' 



(24) 



E 



9nfn 



(25) 



where 3ff(a„) is the real part of a„, c is the Fourier symbol 
of the output operator corresponding to the performance index 
under consideration, and a, g and / are the Fourier symbols 
of the consensus operator Q, and the position and velocity 
feedback operators ( 12 1 respectively. 



These expressions can then be worked out for the variety 
of output operators C representing the different performance 
measures defined earlier. The next result presents a summary 
of those calculations for the six different cases. 

Corollary 3.2: The following are performance measures 
(PI), (P2) and (P3) expressed in terms of the Fourier symbols 
g, f and a, of the operators G, F, and defining vehicular 
formations and consensus algorithms which satisfy assump- 
tions (A1)-(A5). The array O is that of the standard consensus 
algorithm 

1) Consensus 

a) Local Error: 



yloc ^ \ " On 

M/3 ^ man) 



b) Long Range Deviation: 



ylrd ^ _2 



E 



riiH hrid odd, neZf 

c) Deviation from Average: 

ydav 



3fi(a„) 



(26) 



(27) 



I ^ Wn)- ^^^^ 



"Note that in d dimensions, the transformed state vector is of dimension 
2d for each wavenumber n. 



2) Vehicular Formations 
a) Local Error: 

yloc 



1 1 On 

^ n^O, neZ% 



b) Long Range Deviation: 



Vlr'^ = 2d 



E 



m+'-'+Kd odd, nei.% 
c) Deviation from Average: 

ydav _ ^ 1 

2 ^ a f ' 

n^O, n&% i'"''" 



(29) 



(30) 



(31) 



IV. Upper bounds using standard algorithms 

In this section we derive asymptotic upper bounds for 
all three performance measures of both the consensus and 
vehicular problems. These bounds are derived by exhibiting 
simple feedback laws similar to the one in the standard 
consensus algorithm (|2|. In the case of vehicular formations, 
we make a distinction between the cases of relative versus 
absolute position and velocity feedbacks, and derive bounds 
for all four possible combinations of such feedbacks. 

The behavior of the asymptotic bounds has an important 
dependence on the underlying spatial dimension d. For the 
purpose of cross comparison, all of the upper bounds derived 
in this section are summarized in Table U 

For later reference, we note that the Fourier transform of 
the array O in Eq. (j2]) is a quantity that occurs often, and can 
be easily calculated as 



-2dii 



d 

E 



-2/3 



d 

E 



1 — COS 



27r 

7 

TV 



(32) 



A. Upper bounds in the consensus case 

We consider the standard consensus algorithm (fill. In this 
case the array a is exactly O, and thus expression ( [26] l for the 
local error immediately simplifies to 

yloc 1 - 1 



E 



1 = 



which then implies the following upper bound for the individ- 
ual local error at each site 

yloc ^ 

< — . 

M - 4d/3 

Thus, the individual local error measure for the standard 
consensus algorithm is bounded from above for any network 
size in any dimension d. 

The derivation of the macroscopic error upper bounds are 
a little more involved. First we observe that Vj:^'^ < AVf''"''". 



This is easily seen since first, the sums in (27 1 and (28 1 involve 
terms that are all of the same sign (since a„ < 0), and second. 
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that the sum in ( |27) l is taken over a subset of the terms in ( |28) . 
It therefore suffices to derive the upper bounds for V/"". 

We begin with a simphfying observation. Because the arrays 
a we consider are real, their Fourier symbols a have even 
symmetry about all the mid axes of Z^. More precisely 



^(ni,...,nr.,...,nd} 



for any of the dimension indices r. Assume for simplicity 
that N is odd, and define N := {N + l)/2. The even 
symmetry property implies that the discrete hyper-cube 
can be divided into 2'^ hyper-cubes, each of the size of Z^, 
and over which the values of a can be generated from its 
values over Z^ by appropriate reflections. Consequently, a 
sum like (|28ll can be reduced to 



1 _ 2^^ ^ 1 



We now calculate an upper bound on the deviation from 



average measure ( 28 i for the Fourier symbol ( 32 1 of the 
standard consensus algorithm 



v; 



dav 



1 

4^ 

2^ 

4/3 



E 



1 



. Eti(l-cos(fnO) 

N 

1 



< 



32/3 



E 



cos {^Ur 



1 



)) 



7d \ 1 



(33) 



where the first equality follows from reflection symmetry. 



and the inequality follows from (49 1, and noting that the 



denominator is made up of d terms of the form 

2 



2tt 

cos I —n, 



> 



27r 



/V2 



where the inequality is valid in the range Ur E [0, {N — 1)]. 

The asymptotics of sums in Eq. ( [33| are presented in Ap- 
pendixjB] Using those expressions, we calculate the individual 
deviation from average measure at each site 



\rda 



< 



1 



32/3 

32/3' 
1 



j^2-d 

r dhiN'-'^i) 

I log(^) 
N d=l 
< Cd^{ log(7V) d = 2 , 
1 > 3 



d^2 
d = 2 



(34) 



where we have used N < N, and Cd is a constant that depends 
on the dimension d, but is independent of N or the algorithm 
parameter (3. We note that the upper bounds have exactly the 
same form when written in terms of the network size M — 



B. Upper bounds for vehicular formations 

To establish upper bounds in this case, we use a feedback 
control law which is similar to ([S). This law can be most 
compactly written in operator notation as 

U = TqX + TqV + CJoX + foV, 

where Tq is the operator of convolution with the array O 
defined in the consensus problem (|2|. Note that in the multi- 
dimensional case, all signals are d-vectors, and thus Tq above 
is our notation for a diagonal operator with Tq in each entry 
of the diagonal. The last two terms represent absolute position 
and velocity error feedbacks respectively. The first two terms 
represent a feedback where each vector component of Uk 
is formed by a law like ([T]i from the corresponding vector 
components of and Vk and all 2d immediate neighbor sites 
in the lattice. 

With the above feedback law, the closed loop system (j?]) 
has the following expressions for the Fourier symbols of G 
and F 



= 9o 



1 Jn fo H" Oji , 



(35) 



where O is the Fourier symbol ( 32 1. We impose the additional 
conditions that < and /o < since otherwise the 
closed loop system will have an increasing number of strictly 
unstable modes as N increases. When go ^ (resp. fo ^ 0) 
we refer to that feedback as using absolute position (resp. 
velocity) feedback. There are four possible combinations of 
such feedback scenarios. 

We now use these expressions for the symbols g and / to 
calculate upper bounds on performance measures (PI), (P2) 
and (P3) for all four feedback scenarios. We begin with the 
local error ([29]l which in this case is given by 



V, 



loc 



-1 

4^ 



E 



On 



. (5o + 0n)(/o + 0„) 



(36) 



In the case of relative position and velocity error feedback, 
which corresponds to go = and fo = 0, the sum in 
Eq. (36 1 becomes — ^l/0„. This has the same form as 
ydav £q I for the standard consensus problem, and 
thus will grow asymptotically as derived in Eq. ( [34| . For this 
scenario, the final answer is listed as V^°'^ in Table |l] after 
multiplying by the extra 1//3 factor. In the case of relative 
position and absolute velocity feedback, the sum in Eq. ( [36] ) 
becomes J2 ~l/(/o + 0„). Each term is bounded from above 
by -l/(/o + On) < -l/fo since /„ < and 0„ < 0. 
Thus the entire sum has an upper bound that scales like M, 
which yields a constant bound for the individual local error 
once divided by the network size M. An exactly symmetric 
argument applies to the case of absolute position but relative 
velocity feedback. Finally, in the case of both absolute position 
and velocity feedback /o < and go < implying a uniform 
bound on each term in the sum. Similarly the entire sum scales 
like AI and thus is uniformly bounded upon division by the 
network size. All of these four cases for the local error scalings 
are summarized in Table U 
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TABLE I 

Summery of asymptotic scalings of upper bounds in terms of tlie tlie total network size M and tlie spatial dimensions d. Performance measures are 
classified as either microscopic (local error), or macroscopic (deviation from average or long range deviation). There are four possible feedback strategies in 
vehicular formations depending on which combination of relative or absolute position or velocity error feedback is used. Quantities listed are up to a 

multiplicative factor that is independent of M or algorithm parameter f}. 





Microscopic 


Macroscopic 


Consensus 


1/0 


H 


M d=l 
log(M) d = 2 
1 d > 3 


Vehicular Formations 

Feedback type: abs. pos. & abs. vel. 




1 


Vehicular Formations 

Feedback type: rel. pos. & abs. vel. 

or abs. pos. & rel. vel. 




H 


M d=l 
log(M) d = 2 
1 d>3 


Vehicular Formations 

Feedback type: rel. pos. & rel. vel. 


( M d=l 
^{ log(A/) d = 2 
[ 1 d>3 




f A.p d=l 
\ M d = 2 
Ml/3 d = 3 
log(M) d = 4 
[ 1 d > 5 



We now consider the case of the deviation from average 
measure pT| ) which for our specific algorithm is 



^dav 



E 



1 

(3o + 0„)(/o + 6„)' 



When go < and fo < 0, each term in the sum is bounded 
and the entire sum scales as M. Thus, the individual deviation 
from average at each site is bounded in this case. When either 
/o = or 5o = 0, then the sums scale like — 1/ (since 
the other factor in the fraction is uniformly bounded), i.e. like 



the deviation from average in the consensus case (34i 



The only case that requires further examination is that of 
relative position and relative velocity feedback (go — fo — 0). 
In this case 



v 



d 

2 



< 



E 



1 



1 



28 



E 



where the inequality is derived by the same argument used in 



deriving the inequality (33i. Dividing this expression by the 
network size N'^ and using the asymptotic expressions (52 1 
yields 



dav 



(37) 



log(iV) 



where Cd is a constant depending on the dimension d but 
independent of N or the algorithm parameter /3. Rewriting 
these bounds in terms of the total network size M = N'^ 
gives the corresponding entries in Table |I] where the other 
cases are also summarized. 

We finally point out that yj'"'* < due to an argument 

identical to that employed in the consensus case. We thus 
conclude that the upper bounds just derived apply to the case 
of the long range deviation measure as well. 



The role of viscous friction: It is interesting to observe 
that in vehicular models with viscous friction (|6]), a certain 
amount of absolute velocity feedback is inherently present in 
the dynamics. The model (j6]l with a feedback control of the 
form (j8]l has the following Fourier symbol for the velocity 
feedback operator F 

fn = - + fo + On- 

We conclude that even in cases of only relative velocity error 
feedback (i.e. when = 0), the viscous friction term /i > 
provides some amount of absolute velocity error feedback. 
Thus, in an environment which has viscous damping, perfor- 
mance in vehicle formation problems scale in a similar manner 
to consensus problems. These comments are also applicable 
to the lower bounds developed in the next section. 

The role of control effort: A common feature of all the 
asymptotic upper bounds of the standard algorithms just 
presented is their dependence on the parameter /3. If this 
parameter is fixed in advance based on design considerations, 
then the algorithm's performance will scale as shown in 
Table |l] However, it is possible to consider the redesign of 
the algorithms as the network size increases. For example, 
it is possible to increase (3 proportionally to M in consensus 
algorithms to achieve bounded macroscopic errors even for one 
dimensional networks. As can be seen from ([T]), this has the 
effect of increasing the control feedback gains unboundedly 
(in M), which would clearly be unacceptable in any realistic 
control problem. Thus, any consideration of the fundamental 
limits of performance of more general algorithms must account 
for some notion of control effort, and we turn to this issue in 
the next section. 

V. Lower bounds 

A natural question arises as to whether one can design 
feedback controls with better asymptotic performance than 
the standard algorithms presented in the previous section. 
In this section we analyze the performance of any linear 
static state feedback control algorithm satisfying the structural 
assumptions (A1)-(A5), and subject to a constraint on control 
effort. A standard measure of control effort in stochastic 
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settings is the steady state variance of the control signal at 
each site 



^ {K^k} , 



(38) 



which is independent of k due to the spatial invariance assump- 
tion. We constrain this quantity and derive lower bounds on the 
performance of any algorithm that respects this constraint. The 
basic conclusion is that lower bounds on performance scale 
like the upper bounds listed in Table [l] with the control effort 
replacing the parameter /S. In other words, any algorithm with 
control effort constraints will not do better asymptotically than 
the standard algorithms of section |IV-B| This is somewhat 
surprising given the extra degrees of freedom possible through 
feedback control design, and it perhaps implies that it is 
primarily the network topology and the structural constraints, 
rather than the selection of the algorithm's parameters that 
determine these fundamental limitations. 

We now turn to the calculation of lower bounds on both mi- 
croscopic and macroscopic performance measures. For brevity, 
we include only the calculations for the deviation from average 
macroscopic measures. These calculations are a little more 
involved than those for the upper bounds since they need 
to be valid for an entire class of feedback gains. However, 
the basic ideas of utilizing norms are similar, and this 
is what we do in the sequel. In addition, a new ingredient 
appears where the control effort bound, combined with the 
locality property, implies a uniform bound on the entries of 
the feedback arrays. This is stated precisely in the next lemma 
whose proof is found in the Appendix. These bounds then 
finally impose lower bounds on the performance of control- 
constrained local algorithms. 

Lemma 5.1: Consider general consensus Q and vehicular 
formation (j7]i algorithms where the feedback arrays a, g and / 
posses the locality property (A3). The following bounds hold 

Halloo < Ba£{ul) 

llslloo < B, {£{ul]f (39) 

ll/lloo < Bf£{ul), 

where Ba, Bg and Bf are constants independent of the network 
size. 

A. Lower bounds for consensus algorithms 

We start with the deviation from average measure for a 
stable consensus algorithm subject to a constraint of bounded 
control variance at each site 

,2 ' 



£{ul] < W. 



(40) 



We first observe a bound on 5R(a„) that can be established 
from the definition of the Fourier transform 



3?(a„) = ^ 



ttk COS 



27r 



■ k 



E 



ke 

E 

ke 



27r 



-n ■ k 



N 



{-ak) 1 - cos 



2tt 



■k 



where the last equality is a consequence of the condition 
J2keZ'' '^k = 0. For lower bounds on ^l/3fi(— a„), upper 
bounds on 3fi(— a„) are needed. Observe that 



Ofe I 1 — cos 



27r 

TV' 



-n ■ k 



< ^ Iflfcl f 1 - cos 



2ti 



■ k 



< 



k£l 



where the second inequality follows from (48 1. The last 
quantity can be further bounded by recalling the locality 
property ( [TT| i, which has the consequence 

\ak\ {hni H h fcrfTid)^ 

kez% 

= E '"''I ihrii + ■ 

kez%, \k,\<q 

< E l"*^' (^"1 + 

0^keZ%, \k,\<q 

= q'^{ni + ---+ndf ^ |afc|. 

o^feez^, \k,\<q 

Now the locality property can be used again to bound the 



kdnd)^ 



above sum using the the control effort bounds p9[ l and ( 40 1 

< (2g)1|a|U < {2qfBaW. (41) 

O^keZf,., \k,\<q 

Putting the above together gives 

1 



E 

Ar2 



-3?(a„) 



1 



- T:^(2q¥+2 BaW E + + 

> ^ V ^ 



where the last inequality follows from (50i, and Cd is a 
constant independent of N. 



Finally, utilizing ( 5 1 1 and dividing by the network size M = 
N"^, a lower bound on the deviation from average is obtained 



dav 



> 



Cd 

w 



j^2^d 



1 



dhi^-N'-') 



> Cd 



\og{N) 

N d=l 

log(7V) d = 2 , 

1 d>3 



d^2 
d = 2 



(42) 



N 



where by a slight abuse of notation, we use Cd to denote 
different constants in the expressions above. We observe how 
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the lower bounds ( |42] l have the same asymptotic form as the 
upper bounds for the standard consensus algorithm ( [34) l, but 
with the control effort bound W replacing the parameter /3. 

B. Lower bounds for vehicular formations 

We recall the development of the upper bounds for vehicular 



formations in Section IV-B The Fourier symbols of general 



feedback gains G and F have a similar form to ( 35 1, and can 
be written as 



9o 



In, 



In /o 



(43) 



where go, fo and 7, are the absolute and relative feedback 
terms respectively. As before, we impose the conditions that 
9o, fo < 0. We assume that we have a control effort constraint 
of the form (|40]l. 

The case of absolute position and absolute velocity feedback 
has upper bounds which are finite, and the question of lower 
bounds is moot. For the other three cases, lower bounds 



on ( 3 1 1 are established using upper bounds on the symbols 
g and / which can be derived as follows 



ll/lloo < ll/lli < (2g + l) ll/lloo < {2q + l)BfW, 
where the inequalities follow from ([46|, the locality property. 



and ( 39 1 respectively. For g we similarly have 

Consider now the case of relative position and absolute 
velocity feedback. A lower bound is established by 



V 

' 1 



dav 



> 



E 

n^O, nel 

d 



1 



1 



2{2q+\)Bf W E 



\9n 



Now a lower bound on the sum can be established in exactly 
the same manner as (|42]) in the consensus case since g is 
a symbol of a local relative feedback operator. The case of 
relative velocity and absolute position feedback is similar with 
the exception that the factor of is replaced by 

The final case to consider is that of relative position and 
relative velocity feedback. One can repeat the same arguments 



made in the consensus case up to equation (41 1 for both 5„ 
and /„ to state 



Q^n&% \9n\\fn 



> 



ci 



E 



1 



Iffllooll/lloo ^ ^-1^ ("i + --- + nd)4 



1 



where ci and C2 are some constants independent of N and 
W. The asymptotic behavior of this expression (divided by the 



network size) was given earlier in (37 1. We thus conclude that 



the lower bounds in this case are exactly like the upper bounds 
shown in Table [I] for relative position and relative velocity 
feedback, but with the -Jy term replaced by 



1 

W3 • 



VI. Examples and Multiscale Interpretation 

Numerical simulations of cases where macroscopic mea- 
sures grow unboundedly with network size show a particular 
type of motion for the entire formation. In the one dimensional 
case, it can be described as an accordion-like motion in which 
large shape features in the formation fluctuate. Figure [T] shows 
the results of a simulation of a 100 vehicle platoon with both 
relative position and relative velocity error feedbacks. This 
corresponds to a control strategy of the type for which upper 
bounds were calculated in section [TV-B | with go — fo — 0- 
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Fig. 1. Vehicle position trajectories (relative to vehicle number 1) of a 100 
vehicle formation all of which are subjected to random disturbances. Top 
graph is a "zoomed out" view exhibiting the slow accordion-like motion of 
the entire formation. Bottom graph is a zoomed in view showing that vehicle- 
to-vehicle distances are relatively well regulated. 



An interesting feature of these plots is the phenomenon of 
lack of formation coherence. This is only discernible when 
one "zooms out" to view the entire formation. The length of 
the formation fluctuates stochastically, but with a distinct slow 
temporal and long spatial wavelength signature. In contrast, 
the zoomed-in view in Figure [T] shows a relatively well 
regulated vehicle-to-vehicle spacing. In general, it appears 
that small scale (both temporally and spatially) disturbances 
are well regulated, while large scale disturbances are not. 
An intuitive interpretation of this phenomenon is that local 
feedback strategies are unable to regulate against large scale 
disturbances. 

In this paper we have not directly analyzed the temporal 
and spatial scale dependent disturbance attenuation limits of 
performance. However, it appears that our microscopic and 
macroscopic measures of performance do indeed correspond to 
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small and large scale (both spatially and temporally) motions 
respectively. We next outline a more mathematical argument 
that connects these measures. 

Mode shapes: To appreciate the connection between 
norms and mode shapes in our system, consider first a general 
linear system driven by a white random process 

X ~ Ax + w. 

When A is a normal matrix, it is easy to show (by diagonaliz- 
ing the system with the orthonormal state transformation made 
up of the eigenvectors of A) that the steady state variance of 
the state is 

liin£KW.W} . 

where the sum is taken over all the eigenvalues A,; of A. 
Thus we can say that under white disturbance excitation, the 
amount of energy each mode contains is inversely proportional 
to its distance from the imaginary axis. In other words, slower 
modes are more energetic. Now, all the systems we consider 
in this paper are diagonalizable (or block-diagonlizable) by 
the spatial Discrete Fourier Transform. In addition, for the 
standard algorithms, we have the situation that slow temporal 
modes correspond to long spatial wavelengths. This provides 
an explanation for the observation that the most energetic 
motions are those that are temporally slow and have long 
spatial wavelengths. 

String instability: While string instability is sometimes an 
issue in formation control, the phenomenon we study in this 
paper is distinct from string instability. The example presented 
in this section is that of a formation that does posses string 
stability. For illustration, we repeat the simulation but with 
disturbances acting only on the first vehicle. The resulting 
vehicle trajectories are shown in Figure [2] It is interesting 
to note that temporally high frequency disturbances appear 
to be very well regulated, and do not propagate far into the 
formation, while low temporal frequency disturbances appear 
to propagate deep into the formation. What is not shown in 
the figure is that low frequency disturbances are eventually 
regulated for vehicles far from the first. This is consistent 
with the intuitive notion discussed earlier that local feedback 
is relatively unable to regulate large scale disturbances. 
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Fig. 2. Vehicle position trajectories (relative to leader) of the first few of a 
100 vehicle formation. Only lead vehicle is subjected to random disturbances. 
Vehicle trajectories exhibit regulation against that disturbance, indicating the 
absence of string instability. 



Multi-scale properties of disturbance rejection: An intrigu- 
ing explanation of the above example and our scaling results 
is as follows. The macroscopic error measures capture how 
well the network regulates against large-scale disturbances. 
In large, one dimensional networks, local feedback alone is 
thus unable to regulate against these large-scale disturbances, 
and global feedback is required to achieve this. This seems 
rather intuitive. Perhaps surprisingly, in large networks with 
higher spatial dimensionality, local feedback alone can indeed 
regulate against large-scale disturbances. This follows for 
networks for which the macroscopic error measure is bounded 
irrespective of network size. The "critical dimension" needed 
to achieve this depends on the order of the node dynamics 
as well as the type of feedback strategy as shown in Table |l] 
(e.g. dimension 3 for relative position and absolute velocity 
feedback, and dimension 5 for relative position and velocity 
feedback in cases of vehicular formations). 

VII. Discussion 

A. General networks 

The networks considered in this paper are ones which can 
be built on top of a Torus network. Some concepts, such as 
coherence and microscopic and macroscopic errors are easily 
generalized to arbitrary networks. The correct generalization 
of the concept of spatial dimension however is more subtle. 

For any network of dynamical systems for which a distance 
metric is defined between nodes (e.g. from an imbedding of 
the network in M"), the notion of long range deviation can be 
defined as done in this paper. The calculation of that quantity 
involves system Grammians and may even be written in 
terms of the underlying system matrices for certain structures. 
Thus coherence measures can be calculated numerically for 
such networks. However, more explicit calculations to uncover 
scaling laws as network size increases will clearly require 
more analytical expressions for the system norms in such 
networks. 

To generalize the present results, one would require a 
notion of how to grow the network size while preserving 
certain topological properties such as the spatial dimension. 
Preliminary results on self-similar and fractal networks have 
been obtained | [T4| . The proper notion of spatial dimension to 
capture coherence in general graphs remains a research topic 
at this time. 

B. Distributed estimation and resistive lattices 

The results presented here have a strong resemblance to 
results on performance limitations of distributed estimation 
algorithms based on network topology p3) , |T6| , where 
asymptotic bounds similar to ( [34| first appeared in the controls 
literature (see also 1 17 1 where a consensus problem with noisy 
observations is analyzed yielding performance bounds like 
the consensus upper bounds we have in the present paper). 
In that work, the arguments are based on an analogy with 
effective resistance in resistive lattices and certain imbeddings 
of their graphs in d-dimensional space |18|. It is not clear 
how the resistive analogy can be generalized to cover the case 
of second order dynamics (i.e. vehicular formations), or the 
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lower bounds on more general control laws. We have therefore 
avoided the resistive network analogy in this paper by directly 
using the multi-dimensional Fourier transform and reducing 



all calculations to sums of the form ( 5 1 1 resulting in a self- 
contained argument. 

It is interesting to note that the original arguments for the 
asymptotic behavior in resistive lattices p9| in the physics 
literature are based on approximations of the Green's function 
of the diffusion operator in d-dimensions, for which the 



underlying techniques are approximations like (54 1 



C. Order of local dynamics 

We have attempted to keep the development general enough 
that it is applicable to networked dynamical systems whose 
dynamics are not necessarily those of vehicles in any particular 
physical setting. What we refer to in this paper as consensus 
and vehicular formations problems respectively represent net- 
works where the local dynamics (at each site) are first and 
second order chains of integrators respectively. The dynamical 
models are such that the stochastic disturbance enters into 
the first integrator, and the performance objectives involve 
variances of the outputs of the last integrators at each site. 
One generalization of this set up is where the local dynamics 
is a chain of n integrators. It is then possible to show that 
(by retracing the arguments for the vehicular formations case 



and generalizing (52 1) the cutoff dimension to have bounded 



macroscopic measures with only local relative state feedback 
is 1 + 2n. 

D. LQR designs 

It was observed in f20l that optimal LQR designs for vehic- 
ular platoons suffer from a fundamental problem as the platoon 
size increases to infinity. These optimal feedback laws are 
almost local in a sense described by |21|, where control gains 
decay exponentially as a function of distance. The resulting 



optimal feedbacks |20| suffer from the problem of having 
underdamped slow modes with long spatial wavelengths. Thus, 
the same incoherence phenomenon occurs in these optimal 
LQR designs where the performance objective is composed 
of sums of local relative errors (leading to feedback laws with 
exponentially decaying gains on relative errors). 

E. Measuring performance in large scale systems 

In spatial dimensions where performance scalings are 
bounded, the underlying system eigenvalues still limit towards 
zero, suggesting ultimate instability in the limit as M — > oo. 
However, measures of performance remain bounded in these 
cases. In such cases the locations of internal eigenvalues are 
not a good indication of the system's performance in the limit 
of large networks. 

Take the consensus problem over as an example. The 
"least damped eigenvalue" (other than zero) quantifies the 
convergence time of deviation from average (in the absence 
of stochastic disturbances), and it scales as 



as can be shown by explicit eigenvalue calculations | [T2) , 
113]. If we use this quantity as a measure of performance, 
it indicates that performance becomes arbitrarily bad (in the 
limit of large N) in any spatial dimension d. On the other 
hand, consider the use of a macroscopic error measure like 



the variance of the deviation from average (19 1 in the presence 
of stochastic disturbances. That quantity can be expressed in 
terms of the system eigenvalues as 
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—y- 



(45) 



where the sum is taken over all the system's eigenvalues other 



than zero. Note that this sum is just (28 i rewritten to emphasize 
the contrast with (|44]i. 

The important observation is that (|44]i indicates that as 
network size increases, the system eigenvalues approach the 
stability boundary, indicating an eventual catastrophic loss of 
performance in any spatial dimension d. On the other hand, 
( |45] ) is uniformly bounded in dimensions c? > 3 (as shown 
in (|34|i), thus implying well behaved systems as quantified by 
the macroscopic performance measures. A similar point to the 



above has been recently made 1 22 1 



The least damped eigenvalue is traditionally used as an 
important measure of performance. The examples in this paper 
demonstrate that for large scale systems, it is not a very 
meaningful measure of performance, and that the general 
question of how to measure performance in large scale systems 
is a subtle one. 



F. Detuning/mistuning designs 

It is shown in |23 | that spatially-invariant local controllers 
for platoons have closed loop eigenvalues that approach the 
origin at a rate of 0{^). A "mistuning" design modification 



is proposed [23 1, resulting in spatially- varying local controllers 
where the closed loop eigenvalues approach the origin at 
the better rate of 0{j^). In this paper, we have not used 
the real part of the least damped eigenvalue as a measure 
of performance but rather the variance of certain system 
outputs. This amounts to using an norm as the measure 
of performance. It was shown in 121] that for spatially- 
invariant plants, one can not improve performance with 
spatially-varying controllers. The resulting controllers however 
have exponentially decaying gains rather than completely local 
gains. The problem of designing optimal controllers with 
a prescribed neighborhood of interaction remains an open 
and non-convex one. It is an interesting and open question 
as to whether mistuning designs for the measures we 
use in this paper can yield local controllers with better 
asymptotic performance than spatially-invariant ones. It was 
also shown |23| that a mistuning design can improve H°° 
performance for platoon problems. This shows that there 
is perhaps an important distinction between H°° and 
measures of performance for large scale systems. A point that 
is worthy of further investigation. 
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Appendix 

A. Multi-dimensional Discrete Fourier Transform 

We define the Discrete Fourier Transform for functions {fk} 
over by 



where n ■ k :— niki + • • • + Udkd- The inverse transform is 
given by 



A := ^ E A- 



where M = iV''. An immediate consequence of the definitions 
are the following bounds 



ll/IIco < ll/lll, ll/lloo < ^ll/lll. 



(46) 



Let (5 be the Kronecker delta on Z^. It's transform is the 
array 1, which is the array of all elements equal to 1. The 
transform of 1 is MS. We use the symbol 1 to denote the 
array of all elements equal to jj. 

If Tf denotes the circulant operator of circular convolution 
with /, then the eigenvalues of Tf are just the numbers {/„}, 
and consequently the trace of Tf is given by the sum 



(47) 



B. Bounds and asymptotics of sums 

The following facts are useful in establishing asymptotic 
bounds. 

1) For any a; G M and any y e [— tt, tt] 



1 — cos(x) < x^, 

2 , 

1 - cos{y) > V . 



(48) 
(49) 



2) Given d integers ni, . . ., n^, 

(m + • • • + Udf < (2d + 1) (n^ + ■ ■ • + nl). (50) 
Proof: 



E^ 



E-? + EE^ 



Using niUj < (inax{ni,nj}) < nf + n^, we get the 



bound 



E-0 ^ E-? + EE("?+".^) 
< E"? + 2rfE"?- 



3) In the limit of large N , 

E ' 



E 

nGZl 



■,2^2 





d^2 




(51) 


log(iV) 


d = 2 








(52) 


log(iV) 


d = 4 



(53) 



where f{N) « .9(^) is notation for 

c.9(A^) < f[N) < cgiN), 
for some constants c and c and all > iV for some N. 



Proof: We begin with ( 5 1 1. Upper and lower bounds on this 
sum can be derived by viewing it as upper and lower Rieman 
sums for the integral 

It doc -j^ ■ ■ • (j,X (I ^ 



-1 ^ ^ -^d 

over the region A < r < 1 for the lower bound, and A < r < 
^/d for the upper bound. Here A = and the asymptotic 
behavior is determined by the lower limit on the integrals, so 
both upper and lower bounds behave the same asymptotically. 

Using the grid points xi = niA, . . ., Xd — ridA, and using 
the volume increment A'', we get 

2 J L ^2 "^^1 



■ dxd 



A^ E 



1 



((Am 



(And)2) 



E 



(54) 



Now the integral can be evaluated using hyperspherical coor- 
dinates by 

1 



^dxi ■ ■ ■ dxd 



A Jo 



-isin'*-2( 



l) . . . sin(0d_2) dr d4>i 



Cd I r 



d-3 



dr. 



where Cd is a constant that depends only on the dimension 
d (and can be expressed in terms of the volume of the unit 
sphere in M.'^). Evaluating this integral, using A = and 
equation ( [54| gives the result ( [5T] i. 

The proof of (52 1 is very similar to the above, with the 



exception that one approximates the integral of , 2. V 2^2 
^ instead. The details are omitted for brevity. 



C. Proof of Lemma 3.1 

For the consensus problem, the state equation is Q, and 
thus the Lyapunov equation ( 23 1 becomes 



Pn^n 
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where we have used H ~ C (and the choice of C depends on 
the particular performance measure being considered). Since 
all quantities are scalars, this equation is immediately solved 

(A* 

for pn = 2SR(a ) for n ^ Q. la the case n 0, we look at 
the integral definition (22 1, conclude that Cq = implies that 
Po = 0. Thus, the sum (21 1 is calculated to be (|24]). 



For the vehicular problem, the state equation is (j7]i with the 



] . The Lyapunov equation ( 23 1 



" 


Gn 






Z-a 


+ 




Zn 




/ 


K . 




. z: 






. K 







output equation = [ C 
becomes 

/ 

— C * C,j 


where each of the submatrices is of size dxd. From the above, 
we extract the following matrix equations 



F*Y 



Xn 
Zl 



Y F 



Since we are only interested in the quantity 



BjiPnBr 



[0 /] 



Xn 

Zl 



Zn 
Yn 




0. 



(55) 
(56) 



tr K, 



then only equations (55 i and (56 1 are relevant. The coordinate 



decoupling assumption (A5) on the operators G, F and C 
implies that the matrices F„ and Cn are all diagonal. It 
follows that Zn, Xn and Yn are also diagonal, and the above 
matrix equations are trivial to solve. Z„ has the solution 



' 1 /V* *A 

1 L/jiL'Ti- 



Similarly, equation (56i is solved to yield 
1 



Y 



:{GnFn) ^G^Cn 



for n 7^ 0. For the unstable mode at n = 0, the integral (22i 
can be easily evaluated to yield Zq — Q (since Co = for 
all the performance measures we consider). Adding in the 
assumption that all matrices are diagonal with equal elements, 
we obtain in summary the total iJ^ norm of the vehicle 
formation problem (j7|l is given by 



n^O nel 



9nfn 



where the multiplicative factor of d comes from taking the 
trace of Y,. 



Proof of Corollary \3.2\ 

1 ) Consensus: The local error measure output operator C is 

^ Ohy the identity ( 16 1. 
gives the result for V^°'^. 



given by for which C*C 
Combining this with Lemma 



3.1 



The long range deviation measure has the output operator 



defined in (18 1, which has the Fourier symbol 

_ 1 _ „-«ir(niH hnd) 

l^n ± c , 



from which we conclude that 

{ni + ■ ■ ■ + rid) even 
2 [rii + • • • + nj) odd 



Combining this with Lemma 



3.1 



1 n ^n 



gives the result for V^"^"^. 
The deviation from average output operator is C I — T^, 
or equivalently, the convolution operator Tso_^. The corre- 
sponding Fourier symbol is the Fourier transform of the array 
5° — 1, which gives 

n = 
n 7^ ■ 

Putting this in the general formula ( [24| yields the result for 

^ c 

2) Vehicular formations: The derivations for this case are 
very similar to those for the consensus problem and are 
therefore omitted for brevity. 

Proof of Lemma \5.1\ 

We rewrite the dynamics of the consensus algorithm so that 
u is an output 

X = TaX + w 
U = TaX. 

The present task is then to calculate the norm from w to 
u. Applying formula (24i with Tq as the C operator yields 

IT = (-a„), 



E 



n^O, n& 



nei 



after observing that a„ is real and clq = 0. Furthermore, 
our stability condition requires that for all n, (— a„) > 0, 
implying that the sum above is the ^^-norm of {a„}. Putting 
this together with the bound ( |46| ) gives 

Iklloo < ^llalli = 2£{ul). (58) 

In the vehicular formations case, the dynamics are given 
by (|7]i together with the output equation 



[G F] 



X 
V 



Formula ( 25 i is not applicable here since the output depends 
on all states, but the norm for this system can be calculated 



in a manner similar to the proof of Lemma 3.1 They Lyapunov 
equation in this case becomes 



" 


Gn 


Xn 


Zn 




Xn 


Zn 


/ 


K . 


. Z*n 


Yn . 




. K 


Yn . 





Gn 



G nGn 
F*(l 



I 

Fn 

G* F 
F* F 



from which we extract the matrix equations 

G„Z„ + ZnGn ~ 



GnYn 

F*Y 
^n ^ n 



n n 
Xn 

Zl 



G nGn 

'G* F 

^n^ n 

' F* F 



Only the first and last equation need be solved since we are 
only interested in tr(yri). All of the above are d x d diagonal 
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matrices with equal entries, so we solve the equations in terms 
of a single entry as 

1 . 



Zfi C/^ Zji 

'^fnyn — fn '^^n ^ Vn 

The norm of the system is then 



2^ 

-fn + 



1 



E 



-fn + yffn 
In 



l/ll 



7 



.911 1 



where the last equation follows from the stability conditions 
fn < 0, cjn < 0. This inequality has two consequences after 
observing that X^neZ'' tr(F„) = MS {u^} and using (46 1 



l/ll 



\ig\\, < Im8{uI} 



(59) 
(60) 



The second inequality can be used to bound ||g||oo as follows. 
First 

1 



1 M, MM, • , 1 

h .911 1 > 11.911 1 



/II 



An upper bound on ||/||oo is derived from 

ll/IU < ll/lli < (2-z)1|/||oo, 
where the last inequality follows from the locality assumption 



on /. Combining these last two bounds with ( 60 1 yields 



llslli > 



which when combined with p9\ gives 



1 



(29)" 



M\\g\\ 



\\9\\ 



< 



\\f\\^£{ul} < B, {E{ui]y 
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